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Abstract 

Ab initio calculations indicate that topological-defect networks in graphene display the full vari- 
ety of single-particle electronic structures, including Dirac-fermion null-gap semiconductors, as well 
as metallic and semiconducting systems of very low formation energies with respect to a pristine 
graphene sheet. Corrugation induced by the topological defects further reduces the energy and 
tends to reduce the density of states at the Fermi level, to widen the gaps, or even to lead to gap 
opening in some cases where the parent planar geometry is metallic. 
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Since its recent experimental discovery^, graphene has emerged as an unique condensed- 
matter system, for its remarkable electronic and electromechanical properties, as well as a 
test-bed for relativistic-dynamics phenomena, due to the linear and isotropic dispersion of 
the electronic states near the Fermi level, that leads to the massless-Dirac-fermion nature of 
its low-energy electronic excitations^. Devising gap-opening and gap-control mechanisms is 
one of the key issues related to the prospective application of graphene as the new holy-grail 
of electronic devices^. Several mechanisms for changing, perhaps tailoring, the electronic 
structure of graphene have been proposed. Experimental and theoretical works have indi- 
cated the occurrence of gap opening at the Dirac point in graphene due to interaction with 
insulating substrates^. The generation of additional Dirac cones of massless fermions has 
also been identified, due either to imposed periodic perturbations^, or associated with a pe- 
riodic pentagon-heptagon superstructure on grain boundaries in polycrystalline graphene^. 
Of no smaller relevance are the effects of corrugation on the thermodynamic stability and 
transport properties of substrate-supported and free-standing graphene^ 1 ^— . What causes 
corrugation in graphene is an open issue: the possibility of corrugation being an intrinsic fea- 
ture of the material, due to the interaction with a substrate^, or the effect of adsorbates^, 
are scenarios that have been put forward in the literature. 

In this work, we employ ab initio calculations to examine the possibility of topological 
curvature, its stability and energetics, and its effects on the electronic structure of graphene 
sheets. Topological curvature fields are generated when, following the formation of a Stone- 
Wales defect (SW) by a 90° rotation of a carbon-carbon bone- 2 -, the two pentagons and 
two heptagons comprising the SW are separated, and the defect is fully dissociated. We 
consider periodic lattices containing two pentagons and two heptagons, such that these 
isolated odd-membered rings develop local curvature fields, while preserving the null average 
curvature of the parent graphene sheet. Our study broadens the scope of the pentaheptite 
allotrope of carbon previously introduced by Crespi and collaborators^ 2 , that consisted of 
a two-dimensional (2D) carbon sheet composed entirely of pentagons and heptagons. Our 
results raise the possibility of the generation of semiconducting and metallic graphene sheets 
by the introduction of topological defects (TD). We find that periodic networks of TDs 
may present the full variety of electronic structures, at the single-particle level, including 
genuine metals, "graphene- like" null-gap semiconductors, and also finite-gap (~0. 02-0. 75 eV) 
semiconductors. A remarkable result of our calculations is the identification of metallic and 
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semiconducting structures with very small formation energies (~0.1 eV/atom with respect 
to graphene), in which the corrugation induced at the TD sites plays a key role in reducing 
the energy. Being a full first principles address of fully dissociated TDs, our work differs 
from previous theoretical works on corrugated graphene and on topological defects in this 
materiat^M In particular, our results cast doubt on the assumption that Dirac points are 
robust against the introduction of a distribution of TDs in graphene 14 . 

Our calculations are performed in the framework of Kohn-Sham density functional the- 
ory (DFT), within the generalized-gradient approximation (GGA) 1 ^ and norm-conserving 
pseudopotentials in the Kleinman-Bylander factorized form 1 ^. We use the LCAO method 
implemented in the SIESTA code^-, with a double-zeta basis set plus polarization orbitals. 
Calculational parameters are set such that total-energy differences are converged to better 
than 1 meV/atom, with residual atomic forces no larger than 0.02 eV/A. Single-particle 
energy eigenvalues are converged to ~1 meV. 

The pentaheptite form of graphene investigated by Crespi et alr^- is shown in Fig. Ufa). 
The structures we investigate in this work, which can be considered a generalization of the 
pentaheptite, are generated from three different seed structures, as shown in Fig. [TJ The 
first seed in Fig. Ufa) is the pentaheptite itself, with eight atoms in the unit cell. The second 
seed, in Fig. [TJ[b), has a 16-atom unit cell with a SW surrounded by hexagons. The third 
seed, shown in Fig. Hfc), has a 24-atom unit cell where the pentagons and heptagons are 
already dissociated, each being fully surrounded by one shell of hexagons. Larger unit cells 
can be constructed from the seeds, by introducing shells of hexagons surrounding the TDs. 
In order to establish a nomenclature, let the 8-atom seed be called Sn, with a TD number 
density n(l, 1) = 0.18 (in units of A~ 2 ). The 16-atom seed is then labeled S 2 i, and the 
24-atom seed is S31 . From these we can generally build an Sy geometry which is the j-th 
member of the z-th family, having an 8 x i x j 2 -atom unit cell and a TD number density of 
n(i,j) = n(l, X j 2 ). Here, we consider a total of nine planar geometries, of which: four 
are generated from the 8-atom seed, namely Sn, Si 2, S13, and S14; three are built from the 
16-atom seed, namely S21, S22, and S23; and two from the 24-atom seed, namely S31 and 

C 18 
>->32 — • 

For all but the Sn and S 2 i geometries, the SW defects are fully dissociated, and cor- 
rugated structures can be generated by allowing for the local curvature fields generated by 
the isolated TDs to develop. Thus, a total of 16 structures were considered, nine of which 
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are planar and seven are corrugated. Figure [2] shows a 2x2 unit cell of the planar and the 
corrugated S23 geometries. In all cases, the planar geometries did not develop corrugation 
spontaneously under structural relaxation. In order to investigate whether this was an in- 
dication of metastability of the planar forms, we built slightly corrugated initial geometries 
with random off-plane shifts of the atoms imposed on the planar structures. Upon struc- 
tural optimization the S 3l [n(3, 1) = 0.060], S 12 [n(l,2) = 0.045], and S 22 [n(2,2) = 0.023] 
geometries relaxed back to the planar form for off-plane shifts of 0.01 A, and developed 
corrugation for larger shifts, with the exception of the 531, with a higher TD concentration, 
that only developed corrugation for shifts of 0.1 A or larger. Thus, we find numerical evi- 
dence for a metastability of the planar structures with high TD concentrations, down to a 
critical density of TDs (n c ~ 0.023), below which a corrugation instability sets in, as indi- 
cated in Fig. [3j The instability is related to the fact that, for the more diluted structures, 
large incipient conical and saddle-like regions are associated with the isolated pentagons 
and heptagons, respectively, and the cost of flattening such regions onto the plane becomes 
correspondingly larger—. 

In Fig. [3] we show, for the 16 structures, the formation energy (with respect to pristine 
graphene) Ef = E tot /N — Egraph (in eV/atom) as a function of the TD number density, where 
N is the number of atoms in the unit cell, E tot is the the total energy per unit cell, and E graph 
is the total energy per atom for a pristine graphene sheet. Since the planar sheets remained 
numerically stable, we include in the figure the energy of their "relaxed" geometries. Clearly, 
dissociation of TDs, when constrained to the plane, incurs an energy cost that behaves 
linearly with n (with a negative slope), with the exception of the S14, a structure in the 
limit of the connectivity of such geometries in the plane^. Corrugation leads to a reduction 
of the formation energy that is comparable to or larger than the formation energy itself for 
all but the S31 and Si 2 geometries. We define a corrugation energy by the difference between 
the formation energies of the planar and corrugated geometries E c = E^ ane — Ej orr . Ef and 
E c values are included in Table I. Note the correlation between Ef of the planar geometries 
and E c in Fig. [3j the more diluted structures have large Ef values in the planar form and 
also large values of E c , leading to a trend, for the corrugated geometries in each family, of 
decreasing formation energies as the TD concentration decreases. 

Of note also is the fact that the S 2 j family is a low energy branch in Fig. [31 both for planar 
and corrugated forms, while the Sij and S^j families have comparable energies. Among the 
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seed geometries, the S21 has a much lower formation energy than the Su and S31. This 
can be understood by comparing the seed structures: the Su seed has zero hexagons and 
four TDs in the unit cell; the S21 has four hexagons and four TDs, and the the S31 has 
eight hexagons and four TDs. In the latter, the hexagons are highly strained due to the 
SW dissociation, and this structure has the highest Ef among the seeds. The S 2 i has four 
slightly strained hexagons in its periodic unit, and this considerably lowers its energy when 
compared with the Su and S31. As a result of the balance between hexagons and TDs in the 
seed, the bond-stretching and bond-bending distortions in the S2j family are smaller than in 
the other two families. To quantify these structural features, we compute the Keating-model 
elastic energies of the 16 structures. The results are included in Fig. [3] and in Table I. Note 
that the ab initio Ef and the Keating energies follow the same overall trends, with the S2j 
family appearing as a low-elastic-energy branch in Fig. [3j 

Also related to these structural features of the seeds is the degree of corrugation of each 
structure. We define the atomic corrugation Z c of the non-planar structures as the distance 
of the center of the protruding pentagons to the average plane. The values of Z c are given 
in Table I, and the inset in Fig. |3j shows the behavior of Z c as a function of n. Note that the 
S2j family represents a low corrugation branch, due to its comparably less strained planar 
geometry than the corresponding ones in the Sy and S^j families. This also explains why the 
planar S22 with n(2, 2) = 0.023 is metastable, while the the planar S13 with n(l,3) = 0.020 
is unstable against corrugation. 

The diverse nature of the electronic states of these TD periodic lattices can be anticipated 
by examining the DOS plots of the seed structures in Fig. CD The DOS for the Su seed, in 
Fig. QJd), displays a metallic character, with an ab initio value for the DOS at the Fermi 
level that is half of the tight-binding result in Ref. |l2|. The S21 seed, in Fig. [H(e), shows 
a "graphene-like" Dirac-fermion DOS that derives from linear electronic bands crossing 
each other at the Fermi point, forming a Dirac cone on a much reduced energy scale (~[- 
0.1,0.1] eV), beyond which other bands intervene and break the electron- hole symmetry. 
The S31 seed is a semiconductor [DOS is shown in (f)] with a GGA-DFT gap of ~0.6 eV. 

Regarding the general aspects of the electronic structure of all planar and corrugated 
structures in our study, a first observation concerns the role played by topology, relaxation, 
and corrugation, in determining the nature of the electronic states. In order to address the 
effects solely of the topological transformation (TT), we compute the DOS and the electronic 
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TABLE I: Physical properties of TD networks in graphene. Formation energy Et, corrugation 
energy E c (see text), and elastic energies E e i_ are in units of eV/A. Z c (in A) is the measure of 
corrugation (see text). E g (in eV) is the energy gap for the semiconducting systems. The value of 
DOS at the Fermi level (in number of states per energy per atom) is given for the metallic systems. 
The accumulated charge on the pentagons Qp is given for the corrugated geometries. 
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bands of the (planar) seed structures in three different situations, as also shown in Fig. [Q 
(1) a structure with only the TT imposed on graphene, without any relaxation of either 
atomic positions or lattice vectors; (2) with the TT and relaxation of atomic positions only; 
(3) and finally, with the TT and also relaxation of the lattice vectors and atomic positions. 
For the Su, in Fig. QJd), relaxation does not bring any significant changes, and the system 
is metallic in the three cases, with little change in the position of the Fermi level. Things are 
quite different for the S21, m Fig- UX e )> which is a metallic system in the first two scenarios 
above, but develops a Dirac-fermion-like band structure when lattice relaxation is allowed, 
with the Fermi level lying on the Dirac point. For the S31 we only computed the electronic 
states for the cases (1) and (3). The DOS's are shown in Fig. [1(f), where case (1) is seen 
to behave as a metal, with the Fermi level lying at the bottom part of a set of bands that 
start above an energy gap of about 0.5 eV. After full relaxation, enough spectral weight is 
shifted to lower energies for the Fermi level to lye on the gap, and the system to become a 
semiconductor. Hence, the TT by itself does not explain the behavior of the electronic states 
we observe in our calculations, with relaxation effects playing a crucial role. The realistic 
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treatment of the interatomic interactions, provided by the ab initio method, is thus a key 
ingredient. 

In Table I we include the values of the gap for all semiconducting TD networks, and the 
value of the DOS at the Fermi level for metallic ones. Generally, corrugation tends to reduce 
the DOS at the Fermi level, to widen the gap, or even lead to gap opening in some cases 
where the parent planar geometry is metallic. There is one exception, however, in the case 
of the 5*i4 which displays a small gap (0.02 eV) in the planar form and is a metal in the 
corrugated geometry. Overall, a total of nine structures are semiconductors, seven of which 
with gaps ranging from 0.39 to 0.74 eV. Five of the metastable systems in our study have 
formation energies below 0.2 eV/atom (Ef = 0.39 eV/atom for a fullerene molecule with the 
same methodology). Among these, the planar S21 is a Dirac-fermion null-gap semiconductor, 
the corrugated Su is a metal, and the corrugated S22, S23, and S32 are semiconductors. The 
corrugated S23 TD network stands out as a semiconductor with very small formation energy 
(0.09 eV/atom) for a carbon system. 

The low formation energies of the corrugated TD networks in our study, as well as the Z c 
values which are in the range of the corrugation observed experimentally^ 1 ^—, suggest that 
consideration should be given to these structures as a model for the rippled graphene samples 
in the experiments. A feature of corrugated graphene that has been experimentally observed 
is the formation of electron-hole puddles, with a concentration of electrons in the higher parts 
of the samples^. In our corrugated geometries we find that the pentagons at the higher 
parts of the ripples are electron-rich (except for the corrugated S31) as shown in Table I. 
Moreover, the low formation energies indicate that the synthesis of such structures may be 
feasible. We anticipate that kinetics would be the only potentially limiting factor, since 
besides the low formation energies, the barriers protecting such structures from reverting 



back to pristine graphene should be of the order of the 7 eV estimate for the Su in Ref. 



12 



Regarding a pressure-driven transformation of graphite into a three-dimensional bulk of TD 



networks, as envisaged in Ref. Il2|, we note that the areal density of all corrugated structures 
in our study is larger than that of graphene, e.g., one carbon atom per 2.35 A 2 for the S23 
compared with one carbon atom per 2.69 A 2 for graphene, which would favor a pressure- 
driven transformation. We add that very recent experiments have indicated the formation 
of both clustered and isolated pentagon-heptagon pairs in graphene samples obtained by 
reduction of graphene-oxide 2 ^, and corrugation was observed in the areas of the samples 
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containing the defects. Further refinement of such reduction techniques may lead to a viable 
synthesis routes for TD networks. Moreover, this also suggests that other chemistry-based 
synthesis techniques may also be possible 12 . 

In summary, ab initio calculations indicate that the electronic structure of topological- 
defect networks in graphene range from Dirac-fermion null-gap semiconductors to finite- 
gap semiconductors and metals, some of them with very low formation energies, offering a 
possible route for tailoring the electronic nature of graphene sheets. Corrugation induced 
by the topological defects further reduces the energy and modifies the electronic structure 
of the sheets. 
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FIG. 3: (Color online) Ab initio formation energy of planar (filled) and corrugated (empty symbols) 
TD networks, as a function of TD density. Symbols with a diagonal stripe filling show the Keating- 
model elastic energies of the planar geometries. Planar structures to the right of the dashed vertical 
line are metastable. Inset shows corrugation height as a function of TD density. Lines are guides 
to the eye. 
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